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Abstract 

All stellar evolution models for nucleosynthesis require an initial isotopic 
abundance set to use as a starting point. Generally, our knowledge of isotopic 
abundances of stars is fairly incomplete except for the Solar System. We present a 
first model for a complete average isotopic decomposition as a function of metal- 
licity. Our model is based on the underlying nuclear astrophysics processes, and 
is fitted to observational data, rather than traditional forward galactic chemi- 
cal evolution modeling which integrates stellar yields beginning from big bang 
nucleosynthesis. We first decompose the isotopic solar abundance pattern into 
contributions from astrophysical sources. Each contribution is then assumed to 
scale as a function of metallicity. The resulting total isotopic abundances are 
summed into elemental abundances and fitted to available halo and disk stellar 
data to constrain the model's free parameter values. This procedure allows us 
to use available elemental observational data to reconstruct and constrain both 
the much needed complete isotopic evolution that is not accessible to current 
observations, and the underlying astrophysical processes. Our model finds a best 
fit for Type la contributing ^ .7 to the solar Fe abundance, and Type la onset 
occurring at [Fe/H] ~ —1.2, in agreement with typical values. 

Subject headings: Galaxy: evolution, Galaxy: abundances, stars: abundances, 
nucleosynthesis, (stars:) supernovae 



1. Introduction 

Yields from stellar simulations depend on the initial isotopic composition of the star. 
For example, the initial composition is important during hydrostatic burning phases for 
neutron capture reactions on initial metals, affecting odd-z nuclei abundances. In massive 
stars the weak s-process yields are constrained by both the initial CNO abundance, which 
is responsible for providing a neutron source, and also the initial Fe abundance which 
supplies the seeds for neutron capture (Pignatari et al. 2010). The initial Fe abundance 
is also important in intermediate and low-mass asymptotic giant branch (AGB) stars for 
seeding the main s-process yields (Kappeler et al. 1989, 2010; Lattanzio and Lugaro 2005). 
The detailed stellar abundances affect the opacity of the star (e.g., the The Opacity Project 
1995), which in turn will affect the structure as well as mass and angular momentum loss, 
which in turn changes the late stellar evolution. Knowing the initial abundances of heavy 
isotopes is crucial for understanding 7-process abundances, which use s- and r-process 
isotopes as seeds (Rauscher et al. 2002; Arnould and Goriely 2003). Since these ^^-isotopes 
are so rare, any difference in the seed abundances propagates to the resulting 7-process 
yields. Finally, some fraction of the initial stellar composition is not processed and will 
return to the interstellar medium (ISM) unchanged, hence the final abundance pattern will 
directly inherit any poorly estimated initial abundances, but in many cases what is made in 
the stars and what was there initially is difficult to disentangle. 

The purpose of galactic chemical evolution (GCE) is to understand how the abundances 
of the elements and their isotopes evolved from the Big Bang to today, and can be used 
for obtaining the isotopic abundances at any metallicity to use as stellar simulation inputs. 
Traditional GCE models typically split a model galaxy into one or more zones that require 
functional forms for infall and star formation rates (Timmes et al. 1995; Chiappini et al. 
1997; Costa et al. 2009; Henry et al. 2010). Other processes such as interstellar medium 



(ISM) mixing, and galaxy mergers are often not addressed, although there have been 
fairly recent efforts to incorporate mergers in a hierarchical model (Tumlinson 2006). The 
model galaxy is then usually evolved by integrating stellar yields over time, hence these 
models require nucleosynthesis yields from stellar simulations as inputs. The difficulty with 
this approach is that in order to provide self-consistent nucleosynthesis yields, the stellar 
simulations need a complete initial set of isotopic abundances. 

Ideally, such a set would come directly from the GCE simulation. This would require 
knowledge of the complete set of stellar yields for the specific abundances of the GCE model 
at each time step and for each zone: the complicating factor is in reality composition is 
not just a function of metallicity, but also a function of environment, i.e., time and space. 
Instead the inputs for the stellar simulations often use scaled solar abundances or the results 
of some other approximation. Furthermore most GCE models usually evolve only a subset 
of the stable isotopes, such as just the iron peak isotopes (Henry et al. 2010) or everything 
from hydrogen up to the iron peak (Timmes et al. 1995; Kobayashi et al. 2006), or only 
elements. An approach different than GCE modeling to describe chemical abundances has 
been done recently by Ting et al. (2012). They perform a principle component analysis 
on elemental data that takes correlations among elemental ratios and roughly delineates 
different processes responsible for these correlations, as well as identify likely sites and 
metallicity regimes congruent with them. These efforts offer verification of the established 
paradigm concerning many astrophysical processes, and may help constrain others whose 
properties remain incompletely known. 

The approach taken here is complementary to traditional GCE methods, however one 
must take care not to confuse our model with a proper GCE model. We construct an 
astrophysical model of all stable isotopes, based on physical principles for the production 
sites and mechanisms. Effectively, we scale isotopic abundances as a function of a chosen 



model parameter. This isotopic model is then mapped to an elemental model by summing 
the isotopic abundances to their respective elements. We then fit the elemental abundances 
against available observational data to obtain numerical values for the free parameters 
of the model. Our completed model, however, gives the average isotopic history of the 
Galaxy, subject to the approximations employed. The benefit of this approach to isotopic 
GCE is that it is not necessary to know or model dynamic and galactic evolution processes 
employed in traditional GCE models such as infall, ISM mixing, and galaxy mergers, whose 
uncertainties are poorly constrained. 

Compared to full GCE calculations, however, our approach is rather simplistic and 
approximative, precisely because we do not integrate stellar yields or address dynamic and 
galactic evolution processes employed in traditional GCE models. The intention here is 
to improve upon the typical standard of scaling isotopic solar abundances by a constant 
factor, which effectively treats all istopic production as primary. Our improvements to this 
standard is to incorporate secondary processes and Type la contributions with separate 
scalings, to better approximate their relative values at desired input metallicities for 
nucleosynthesis studies, andthe resulting isotopic histories from our model can be used as 
inputs in stellar models, and comparison to abundances from other sources like Damped 
Lyman-a systems and dwarf galaxies. This represents an improvement over the previous 
standard of guessing ad-hoc assumptions for the interpolation between the known endpoints 
of solar and big bang nucleosynthesis (BBN), or the use of scaled solar abundances. This 
model could also be seen as a first iteration step for more sophisticated first principle GCE 
models, but at a higher order approximation than just scaled solar abundances, but not a 
replacement for them. Hence while the approximations and comparisons of our model are 
sufficient for this purpose, they would be quite unsatisfactory for describing GCE itself. 

This paper has the following outline. Section 2 introduces the astrophysical processes 
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and sites considered by the model, and the separation of the solar isotopic abundance 
pattern into contributions from these processes is discussed. In Section 3, the model itself 
is introduced and the scaling of the processes as a function of our parametrization is 
explained, using their relative solar contributions and the BBN abundance pattern as fixed 
boundary conditions. This also defines the relevant free parameters for fitting the elemental 
abundances with stellar data. In Section 4, the elemental model is fit to available data and 
best fit parameter values are found. In Section 5, the resulting elemental model is discussed 
and additional results given by the model. The final Section 6 addresses constraints of the 
model and discusses possible extensions for future work. 

2. Astrophysical Processes and Solar Abundance Decomposition 

Our model attempts to cover the essential key astrophysical processes responsible for 
the production of the isotopes in the Galaxy. We first use these processes to decompose 
the solar system abundances. Here the isotopic solar abundance pattern was taken from 
a new data set by Lodders et al. (2009) which gives updated values relative to their 2003 
publication (Lodders 2003). In the following we organize the processes roughly by the mass 
range of isotopes to which they contribute. 

2.1. Big-Bang Nucleosynthesis 

When the Universe was less than ~100 seconds old, protons and neutrons were 
in thermal equilibrium with each other by weak interactions with neutrinos. Upon 
continuing expansion, however, the temperature dropped sufficiently to "freeze out" neutrino 
interactions. At this "freeze-out" temperature, corresponding to kBT~0.8MeV, the weak 
interaction rate became slower than the hubble expansion, and the neutron-to-proton ratio 
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(after some subsequent /3 decay) was fixed to n/p ^ 1/7 (Yao et al. 2006). At this point 
BBN began with deuterium formation, and proceeded to produce non-negligible abundances 
of ^H,^He, ^He, ^Li, and minute abundances of ®Li and isotopes up through oxygen. 

We use the theoretical BBN abundance pattern from Cyburt et al. (2001) and provided 
by Fields (2002). Isotopic contributions from this BBN pattern are taken for ^H, ^H, 
^He, ^He, and ^Li. The remaining negligible contributions of ^Li and the isotopes heavier 
than '''Li are not used. A constant ^Li abundance at low metallicities was observed that 
suggested a primordial abundance, but whose value was much smaller than the predicted 
BBN abundance (Spite and Spite 1982). This "Spite plateau" remains unexplained, and for 
the present model we use the theoretical BBN abundance for ^Li. 

2.2. Light Isotopes 

2.2.1. Helium 

The remaining helium not made during BBN is a product of hydrogen burning and 
is scaled as a primary process. The "release time scale" of helium can vary depending on 
the stellar mass, and relative to metallicity it may scale slightly slower than a true primary 
process (the same is true for C and N). We do not account for this delay in the present 
model. 

2.2.2. V- Process 

The "light" //-process involves interactions among neutrinos and lighter nuclei in CCSNe 
environments (Woosley et al. 1990; Yoshida et al. 2004; Heger et al. 2005). The neutrinos 
elevate nuclei to excited states, which then decay by nucleon emission (Hartmann et al. 



1991). The target nuclei for these interactions that produce hght isotopes are CNO isotopes 
made from hydrogen and hehum, hence this production is primary. This process produces 
^^B (Heger et al. 2003), and some ^Li (Prantzos 2010). The "heavy" //-process also involves 
neutrino interactions, but with target nuclei made from either the s- or r-processes, and is 
responsible for heavy nuclei production such as ^^°Ta and ^'^^La. Due to the requirement 
of pre-existing s- or r-process metals to serve as the target nuclei, the "heavy" //-process 
behaves like the 7-process with respect to metallicity, discussed below in Section 2.5.6. 
Hence we do not distinguish between the "heavy" //-process and 7-process contributions. In 
our model, the //-process stands for the "light" //-process, and the 7-process includes the 
"heavy" //-process. 

2.2.3. Galactic Cosmic Ray Spallation 

Galactic Cosmic Ray (GCR) spallation events occur when energetic protons or 
a-particles impact on existing CNO nuclei in the ISM (Reeves et al. 1970; Meneguzzi et 
al. 1971). GCR spallation contributes to ^Li, ^Be, i°B, and "B (Prantzos 2007, 2010). 
Since spallation occurs on pre-existing CNO nuclei in the ISM, this process is traditionally 
considered secondary. Observations, however, show a primary dependence on metallicity for 
^Be (Prantzos 2007), which is in conflict with the understanding of traditional spallation 
events. Prantzos (2010) proposes a solution to this problem, and posits that GCRs 
accelerated by the winds of rotating massive stars could be abundant in ^Be. If these GCRs 
then hit the forward shock of subsequent SNe explosions (which are proton rich), this would 
satisfy the condition of a primary event. We adopt this proposed solution, and assign 
primary GCR spallation as a mechanism for ^Be production, and secondary GCR spallation 
for producing ^°B, ^^B (Prantzos 2007), and ^Li (Prantzos 2010). Note the only other way 
to explain the primary behavior of ^Be in traditional spallation events is by assuming that 



the CNO content in GCR spallation events is always the same (Duncan et al. 1992). 

2.2.4- Classical Novae 

White dwarfs accreting material from a companion star can undergo outbursts powered 
by thermonuclear runaway in the accreted layer (Truran 2002; Jose and Hernanz 2008). 
Nucleosynthesis occurs on the accreted material that is rich in H and He and dredge-up 
of primary CNO (and ONeMg for ONeMg novae) into the envelope, hence this process 
is primary. Unlike most other primary processes that immediately begin to enrich the 
ISM however, there is likely some onset timescale for novae contributions similar to Type 
la SNe. We do not consider this delay for novae in our model. Simulations have shown 
differing isotopic production below the iron peak, depending on the composition of the core 
and its mass (Gehrz et al. 1998; Jose and Hernanz 2007, 2008). Hence precise abundance 
determinations are difficult to isolate. Many CO novae simulations show production of ^Li, 
13q^ ^^N, ^^O, and ^^F that dominate the ejecta, whereas ONeMg novae additionally show 
contributions to other metals up to ^°K (Jose and Hernanz 2007). We take all contributions 
beyond ^Li to be negligible compared to massive star contributions. This approximation 
holds well for CNO isotopes that have large contributions from massive stars, but for 
isotopes such as ^^F the approximation is less than ideal. In fact it believed that ^^F may 
also be produced in the //-process in core-collapse supernovae (CCSNe; Woosley and Haxton 
1988), during hydrostatic nucleosynthesis in He shell of thermally pulsating asymptotic 
giant branch stars (TP-AGB; Forestini et al. 1992), and in the He core of heavy mass loss 
Wolf-Rayet stars (Meynet and Arnould 2000). The decomposition from these sources is still 
a matter of debate (Abia et al. 2010; and references therein), and we do not address this 
complication in the current model. 
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2.2.5. Light Isotope Decompostion 

Identifying the precise (non-BBN) Li, Be, and B contributions presents a challenge. At 
present, there is no consensus for explaining the solar abundance pattern for the isotopes 
of these elements using the processes that could be responsible. Due to the difficulty in 
determining what fraction each actually contributes to the light isotope solar abundance 
pattern, novae, the //-process, and primary GCR spallation are placed into a single category 
due to their shared primary nature. Traditional GCR spallation is treated separately since 
it is secondary. The relative solar abundance decomposition between these two categories 
for ^Li and ^Li is estimated from Prantzos (2010). The 60% of the solar ^^B is assigned 
to the GCR category and the remaining 40 % is believed to be from a primary source, 
possibly the //-process (Woosley et al. 1990; Prantzos 2007). This additional primary source 
is needed in order to obtain the observed meteoric ^^B/^°B ratio (Prantzos 2010). 

2.3. Low and Intermediate-Mass Stars 

Stellar winds from low and intermediate-mass stars are rich in C and N isotopes 
(Arnould et al. 2003), and provide significant contributions to the solar abundances for 
these isotopes. Fitting the contributions from these sources to data is problematic, since 
contamination from massive stars is always present and difficult to separate out. Hence we 
do not independently address stellar wind contributions in our model and instead combine 
their contributions with those of massive stars, both of which are primary processes. 

2.4. Intermediate-Mass and Iron Group Isotopes 

Hydrostatic burning in massive stars (~10 — 100 M0) synthesizes most isotopes from 
helium up to the iron peak (Burbidge et al. 1957; Rauscher et al. 2002). Stellar winds can 
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eject some of this material over the star's hfe but the explosive stellar death dominates 
the metal yields. CCSNe likely produces between 1/3 and 2/3 of the solar abundance iron 
peak isotopes (Timmes et al. 1995). They also produce the majority of the alpha isotopes 
and many of the intermediate isotopes from ^^0 to the iron group. Type la SNe are 
thermonuclear explosions of accreting white dwarfs (eg., Nomoto et al. 1997; Hillebrandt et 
al. 2000; Woosley 2001), and primarily provides the remaining iron peak solar abundances, 
with some enrichment of other metals (Travaglio et al. 2004b; Maeda et al. 2010). In 
fact, simulations show Type la production of minute trace contributions to the isotopes 
below the iron peak (Nomoto et al. 1997), with the exception of ^"^K. Both CCSNe and 
Type la produce their isotopic yields explosively, which destroys much of the initial metal 
composition. The evolutions of their isotopic products are considered primary. 

Yields for Type la supernovae were taken from the W7 model (Nomoto et al. 1997) 
for isotopes with mass numbers 12 < ^ < 56. The category of "massive star contributions" 
is defined in this context to be the collection of all primary isotopic production with mass 
numbers 12 < yl < 68 not attributed to Type la SNe. This includes all isotopic enrichment 
to the ISM driven by massive star stellar winds and production from CCSNe, the r-process, 
the //-process, novae yields, and stellar winds from low and intermediate-mass stars, with 
production from CCSNe dominating the isotopic abundances in this category. The solar 
contributions from massive stars were taken from the yields of a massive star simulation 
(Heger and Woosley 2010) fitted to a low metallicity data set (Cayrel et al. 2004). The 
simulation included stars in the mass range 10 — 100 Mq, with a Salpeter initial mass 
function (IMF) and a relatively low mixing of 0.1 (Joggerst and Heger 2008) employed 
in a running boxcar method (Heger and Woosley 2010). The energetics of the supernova 
explosion are parametrized by. 
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where E is the explosion energy as a function of the stellar mass M, the explosion 
exponent parameter E^xp (taken as 0.0), and Eq (taken as 1.2 B). This parametrization can 
be motivated by noting that as the mass of the star increases, so does the energy required 
to the lift this mass in the explosion. 

The fit of the yields to the Cayrel et al. (2004) data set gave a x value of 1.802 (see 
Heger and Woosley 2010 for fitting procedure). The heavier massive star contributions 
for mass numbers 57 < ^ < 68 were not taken from Heger and Woosley's massive star 
simulation, and were instead calculated as residuals from the main and weak s-processes, 
discussed below in Section 2.5. 

Under the assumption that Type la are responsible for some fraction / of the observed 
solar ^^Fe abundance, each W7 yield was scaled to this fraction. The scaling factor is given 
as: / ■ X^/X^Q, where X^ is the solar abundance of ^^Fe, and X^q is the W7 yield for ^^Fe. 
Hence each isotopic abundance was scaled by this factor, which shifts the entire abundance 
pattern until the yield for ^^Fe is equal to / ■ X^. The fraction / represents a free parameter 
in the model, which is determined by fitting the elemental scalings against available data 
(in Section 5). The massive star yields were scaled to the remaining contribution to solar 
^^Fe not accounted for by Type la. This factor is given as: (1 — /) ■ X^/X^^^^"^, where 
jj-massivc jg ^j^Q massivc star yield for ^^Fe. An additional scaling of both types of SNe data 
in the range 12 < ^ < 56 was then required to ensure that the massive star and Type la 
contributions summed to the solar abundance for every isotope, as the first scaling using / 
would only guarantee that ^^Fe satisfied this requirement. This additional scaling preserved 
the ratio of each isotopic contribution between the W7 and massive star yields. 
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where Xj q, Xi^^ are the original and (scaled) fitted abundances of isotope z, for either 
the massive or Type la contributions (denoted as superscripts), and X® is the solar 
abundance of isotope i. Note that for clarity this procedure has been explained using two 
successive scalings, but in practice this can be done with a single scaling that achieves both. 
This second scaling preserves the isotopic ratios across each model, but the ratios within 
each model undergo some distortion. That is, the overall abundance patterns of the W7 
and massive star simulations are altered. Nevertheless, the final abundance patterns we use 
still show Type la contributing mostly to the Fe peak, and massive stars contributing to 
CNO and a-elements up to the Fe peak. 



2.5. Heavy Isotopes 

2.5.1. Weak S- Process 

The s-process is one of the four trans-iron processes for making the heavier nuclides 
distinguished here. It synthesizes isotopes via slow neutron capture (relative to the beta 
decay rate; Burbidge et al. 1957). It is responsible for approximately half the heavy isotopes 
beyond iron (Pignatari et al. 2010; Heil et al. 2007). Since this process is characterized by 
neutron capture rates that are slow compared to the beta decay rate of the target nucleus, 
production proceeds along the path of isotopic stability, with ^^Fe playing the role of the 
seed nucleus. In practice many metals could seed this process, but ^^Fe will stand as the 
sole target nucleus due to its large abundance and neutron capture cross section relative to 
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other potential seeds. The weak component of the s-process occurs in massive stars, during 
convective core He burning and shell C burning (Pignatari et al. 2010). The neutron source 
in core He burning is from the reaction ^^Ne(a,n)^^Mg. The ^^Ne nuclei are produced from 
the burning of ^*^N made previously in the CNO cycle. Subsequent C shell burning can 
produce neutrons by itself and provide a-particles to reignite the ^^Ne(a,n)^^Mg neutron 
source. These neutrons are captured on initial ^^Fe present in the star. 

The weak s-process component can only synthesize isotopes along the line of stability 
up to a mass number of A^ 100 (Raiteri et al. 1993), due its smaller neutron exposures 
relative to the main component. Both components of the s-process rely on the existence 
of metals that formed the initial composition of the star, and therefore are considered 
secondary processes. The neutron source for the main s-process, however, is primary (as 
well as some of the shell C burning contributions to the weak s-process), which could result 
in a behavior between that of a primary and secondary process. 

For both the main and the weak components of the s-process, bottlenecks exist at 
closed neutron shells, allowing the abundances at these corresponding mass numbers to 
accumulate into peaks. Three such peaks exist at approximate mass numbers 88, 138, and 
208. Elements of interest at these peaks for representing the scaling of s-process elements 
include strontium, barium, and lead. Good recent compilations of the main s-process can 
be found in Lattanzio and Lugaro (2005), and a good review of the s-process can be found 
in Kappeler et al. (2010). 

Modeling the weak s-process has been less successful than models for the main 
component (Pignatari et al. 2010). Since flow equilibrium is not reached during the 
neutron exposure, uncertainties in the neutron capture cross sections affect the yields of all 
subsequent isotopes (Heil et al. 2007). We need to reproduce the solar system abundances 
for the s-only isotopes for both the main and weak components. For the weak s-process, the 
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s-only isotopes '''^Ge, ^^Se, *°Kr, ^^Kr, ^^Sr, and ^^Sr owe their solar abundances to both the 
main and weak components (Raiteri et al. 1993), and the relative contributions from each 
represent a poorly known a priori constraint. Additionally, recent stellar model calculations 
for the weak s-process have difficulties in both producing sufficient s-only isotopic yields 
without at the same time overproducing many other isotopes beyond their solar abundances 
(Pignatari et al. 2010). 

For the purposes of making a reasonable assessment of the weak s-process contributions 
to the solar system abundances, we performed a calculation using the updated online 
MACs compilation from the KADoNIS project www. kadonis.org (Dillmann et al. 2009). 
The branching points addressed in our calculation that lie along the weak s-process path 
are the unstable nuclei ^'^Cu and ^°Br, with corresponding P'^ and /3~ thermal branching 
ratios calculated from Takahashi et al. (1987). The branching point at ™Se was ignored, 
and was taken via (3~ decay to '^^Br prior to any neutron capture (which would give a ^°Se 
population). This approximation is due to the high temperature dependence of the /3~ 
decay rate of '^^Se (Heil et al. 2008; Makinaga et al. 2009; Waher et al. 1986). ^^Zr was 
treated as stable for this calculation, since the thermal (3^ decay branching ratio is small 
relative to neutron capture rate (Takahashi et al. 1987). After the neutron exposure has 
ended, the remaining ^^Zr abundance is taken to ^^Nb. The ^^Rb branching was omitted 
and it was assumed that it entirely (3~ decays to ^^Sr, again because the neutron capture 
branching ratio is negligible (Takahashi et al. 1987). 

The differential equations for the weak s-process abundances under the classical 
approximation (Kappeler et al. 2010) were solved numerically. We have the following 
system of linear differential equations: 
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where X is a vector of the isotopic abundances, X= (Xfcsg? ^Fcst, • • • , -'^Ruioi)! t is 
the neutron exposure, and A is the matrix of Maxwelhan averaged neutron capture cross 
sections and branching ratios for beta decay, with components Ay = — o", for j = i, with 
Aij_i = +0"j, Aj+ij = (XiPi, and Aj+2,j = (^[{l — Pi), where Pi is the beta decay branching 
ratio for isotope i at branching points (where apphcable). Note (1 — /3j) is the /3^ branching 
ratio, and all other elements in the sparse matrix are zero. 

In solving this system, we did not assume a continuous neutron exposure distribution. 
Instead, different values for single exposures were used. A linear combination of the yields 
from different single neutron exposures was determined that best fits the abundances of the 
s-only isotopes in the weak s-process range. In addition, to ensure that the sum of s-only 
isotopes for the weak and main components were not overproduced with respect to solar, 
the abundances to which the s-only isotopes were fit were taken as the residuals of the main 
s-process yields subtracted from the solar abundances. In addition, 15 % of the solar ^''Kr 
and 3 % of the solar ^^Kr were attributed to the z/]9-process (Kappeler et al. 1989; Pignatari 
et al. 2010) and likewise subtracted. 

We found a linear combination of five neutron exposures that reproduces the desired 
s-only isotopic abundance residuals for ^"Ge, ^°Kr, and ^^Sr Fig. 1. The remaining s-only 
isotopes, ^^Se, ^^Kr, and ^^Sr were under-produced. The only over-productions that resulted 
were for ^^Cu and ^^Y. We found the neutron exposures to be: 0.2, 0.25, 0.3, 0.35, and 
0.4 mb~^. The respective coefficients that the abundances from each of these exposures were 
weighted were: 1.5, 0.2, 0.4, 0.6, and 0.7. We attempted to fit as many s-only isotopes to 
their residuals without generating over-productions, and it was not possible to fit all s-only 
isotopes without causing additional isotopic over-abundances. The over-productions of ^^Cu 
and ®^Y as well as the under-productions of ™Se, ^^Kr, and ^^Sr were scaled with the main 
s-process abundances, to fit with the solar abundances. This scaling was equivalent to the 
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scaling done for massive stars and Type la yields: the isotopic ratios were preserved. 
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Fig. 1. — The calculated weak s-process contributions to the solar abundance pattern. Red 
boxes: s-only isotopes. Black x's: isotopes with contributions from the r- and s-processes. 
The yields for the s-only isotopes ^°Ge, ^°Kr, and ^^Sr are reproduced well, with underpro- 
ductions for ^^Se, ®^Kr, and ®^Sr. The only over-productions were for the non s-only isotopes 
65 Cu and 89 Y. 

In addition to producing isotopes along the path of stability, there is indication that 
the weak s-process component in massive stars is responsible for seeding a non-negligible 
p-isotopic production component. The massive star yields from Rauscher et al. (2002) show 
significant isotopic productions for ^^Se, ^®Kr, and ®^Sr, which cannot be accounted for by 
the z/p-process (since their stellar models include weak s-process and 7-process reactions, 
but do not include the z/p-process). A mechanism for producing these ^^-isotopes are 
7-process reactions on weak s-process seeds, where the weak s-process would enhance the 
abundances of unstable proton-rich isotopes which would then undergo photo-disintegration 



events to give stable ^^-isotopes. The scaling of these ^^-isotopic yields should thus track 
the weak s-process, rather than the traditional 7-process (as far as its dependence on 
metallicity), since the abundances are made in-situ from existing weak s-process seeds. 
It is unclear exactly how much this weak s-process enhanced p-process (WSEP) should 
contribute to the solar ^'^Se, ^®Kr, and ^^Sr abundances. We decided to attribute half of the 
solar abundances to each the WSEP isotopes and the z/j»-process for these three j>-isotopes. 

The uncertainties in our calculation are constrained by the errors of the MACs given 
in the compilation from the KADoNIS project www. kadonis.org (Dillmann et al. 2009), and 
vary by isotope. We do not propagate the uncertainties through our equations, since our 
treatment of the weak s-process is only approximate. 

2.5.2. Lighter Element Primary Process 

Indication for the need of an additional primary process distinct from the r-process 
appears to be implied by ultra- metal poor (UMP) stellar abundances, and was first 
implemented by Qian and Wasserburg (2001) in a two component phenomenological model. 
This process, sometimes referred to as the weak r-process (Truran and Cowan 2000) or 
charged-particle reaction process (Qian and Wasserburg 2007), was named in more general 
terms by Travaglio et al. (2004a) as the lighter element primary process (LEPP), and is 
needed to explain an observed excess of some lower mass {A < 130) elements, notably Sr, Y, 
and Zr, that can not be accounted for by neutron capture processes, photo-disintegration, 
or CCSNe. Investigation of the triple-a and ^^C(Q;,n)^®0 rates indicates that their present 
2o" uncertainty can not account for the necessary production of Sr, Y, and Zr by the weak 
s-process in massive stars (Tur et al. 2009). A more recent nucleosynthesis calculation by 
Arcones and Monies (2011) suggests the interesting possibility of the needed abundances 
being produced in the neutrino-driven winds of ultra-metal poor (UMP) CCSNe, although 
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their yields suffer over-productions of additional isotopes in order to provide the necessary 
Sr, Y, and Zr. In our model we do not separate out the yet unknown LEPP process, but 
instead the abundances produced by this mechanism are absorbed into the massive star 
category. 

2.5.3. Main S-Process 

The main component of the s-process occurs in the thermally pulsating AGB stellar 
phase for stars with M < I.SMq. During hydrogen burning, protons are thought to mix 
downward into the helium layer, which can then be captured on synthesized ^^C to form a 
^'^C pocket (Busso and Gallino 1999). After the subsequent helium flash, a-particles are 
convectively dredged up through this pocket initiating a neutron source via ^^C(a;,n)^^0. 
This neutron source drives the main s-process by capture on pre-existing metals contained 
in the star throughout several helium flash cycles. Due to the longer neutron exposures 
operative in the main s-process (relative to the weak component of the s-process), the main 
isotopic contributions are isotopes with mass numbers ^1 > 88 (Iben 1975; Truran and Iben 
1977). 

The main s-process contributions to the solar abundance pattern were taken from 
Bisterzo et al. (2011), and these yields were re-normalized to the Lodders et al. (2009) 
abundances. Many of the s-isotopes with mass numbers A > 88 were consistent with their 
solar values, however, some were under-produced. For the under-produced s-only isotopes 
that fell within 0.1 dex below their solar abundance, we assumed them to be their nominal 
solar values. This cutoff at 0.1 dex is admittedly somewhat arbitrary. The uncertainties 
in neutron capture cross-sections and modeling of the AGB stars introduces error, yet 
remarkably so many s-only isotopes are reproduced by Bisterzo et al. (2011) to within 25% 
(~ 0.1 dex) of their solar values. Our choice of this cutoff value acknowledges the existence 
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of errors that are often difficult to enumerate, while assuming that since so many s-only 
isotopes are reproduced close to their solar values, the model is likely reliable. 

The resulting abundance pattern is given in Fig. 2. Below a mass number oi A = 88, 
the s-only isotopes show an explicit drop in their production, and the weak s-process has 
to contribute in this region to reproduce the needed solar abundances. For mass numbers 
above A = 88, three s-only isotopes are under-produced (by more than O.f dex) by Bisterzo 
et al. (20f f): ^^^Gd, ^^''Os, and ^^^Pt. We attributed the residual abundances of these 
three isotopes to the 7-process. This treatment of the residuals is very approximate, and 
while ^^^Gd and ^^^Pt lie on the proton-rich side of the stability line, ^^^Os is preceded by 
the more proton rich ^^^Os and ^^^Os, both of which would be more likely candidates for 
7-process enrichment. The contribution of ^^^Gd to elemental Gd is only ~0.2 %, hence 
relatively little 7-process enrichment is needed to populate the residual abundance. The 
isotope ^^^Pt contributes ~0.79% to elemental Pt, however the j>-isotope ^^°Pt contributes 
even less at ~0.01 %, so it is unlikely that a residual abundance from the 7-process would 
populate the needed ~0.7% ^^^Pt without also raising the abundance of the more rare 

190Pt. 
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Fig. 2. — Adopted main s-process abundances from Bisterzo (2011), relative to the solar 
abundance pattern. Red Boxes: s-only isotopes. Black x's: isotopes with contributions 
from the r- and s-processes. Under-produced s-only isotopes that fell within 0.1 dex (dashed 
line) below their solar abundance were assumed to be their nominal solar values. 

The main s-process produces three characteristic peaks at closed neutron shells, and 
we sub-divide the main s-process abundances into Is, hs, and "strong" components. The 
Is component is taken to be all abundances from Fig. 3 up to and including Sr, the hs 
component is taken to be all abundances after Sr and up to and including Ba. The strong 
component is assigned to Pb and Bi isotopes, and is discussed in more detail below. Using 
the solar main s-process abundances given in Fig. 3 as a starting point, we scale the Is, hs, 
and "strong" components separately due to the different neutron exposures required for 
their production. 
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2.5.4- Strong S- Process 

Originally a strong component of the s-process was introduced to address the 
underproduction of Pb not accounted for by the main component at solar metallicities 
(Clayton and Rassbach 1967), and it was believed that a third type of neutron exposure 
was needed to generate the remaining Pb. Sufficient production was later found to be 
present in the low metallicity regime of the main component, and the strong s-process 
component has since been re-interpreted as a low metallicity effect of the main component 
(Gallino et al. 1998). As the metallicity (and hence Fe) decreases the neutron-to-seed ratio 
increases, providing a sufficient neutron exposure to make the heavy Pb and Bi isotopes. 
Simulations of low metallicity AGB stars show production of "strong" s-process isotopes at 
[Fe/H] = —2.6 that exceed production at solar metallicities by several dex, depending on 
^■^C pocket efficiencies (Bisterzo et al. 2010). This implies nearly all of the solar "strong" 
s-process abundances are made at low metallicities. 

We shall continue to refer to Pb and Bi s-process contributions as coming from the 
"strong" component, since they are scaled distinctly from the hs and Is parts of the main 
s-process, but it should be noted this is a convention introduced for clarity, and the strong 
s-process is indeed the low metallicity regime of the main s-process, having a distinctly 
stronger neutron-to-seed ratio. 

2. 5. 5. R-Process 

The r-process synthesizes isotopes beyond the iron peak using rapid neutron capture 
(relative to the beta decay rate; Burbidge et al. 1957). The location of this process has not 
been universally accepted, but historically was first thought to occur in CCSNe environments 
(Hoffman et al. 1997). More recently it has been postulated to occur in shocked surface 
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layers of 0-Mg-Ne proto-neutron stars (Ning et al. 2007; Qian and Wasserburg 2007), and 
also simulations have shown success in reproducing r-process signatures from //-driven 
nucleosynthesis in the He-shell during CCSNe in the low metallicity regime (Z < —3.0; 
Banerjee et al. 2011). A recent principle component analysis has linked r-process elemental 
abundances with alpha-elements, further suggesting CCSNe as a possible site for this 
process (Ting et al. 2012). 

The r-process proceeds far beyond the neutron-rich side of stability and (as with 
the s-process) also bottlenecks at closed neutron shells. Since the closed neutron shells 
are encountered on the neutron rich side of stability, the proton number is lower at these 
bottlenecks for the r-process than it is for the s-process, hence after decay to stability 
the r-process peaks always occur at lower mass numbers relative to the s-process peaks, 
at approximately 80, 130, and 195. Elements of interest for representing the scaling of 
r-process elements that have r-only isotopes include germanium, europium and platinum. 
The explosive environment of the r-process does not directly depend on the initial metallicity 
of the star, and therefore is primary (stellar populations that could be responsible for the 
r-process could depend on Z). 

The r-process solar abundance contributions were determined using the residual 
method, in the mass range 69 < ^ < 238. The r-process contributions in the range 
56 < ^ < 68 were not differentiated from the "massive star category." This choice was 
made to facilitate the decomposition of the solar abundance pattern in this range, where 
both CCSNe yields (from a-rich freeze-out, the a-process, and previous nuclear burning) 
and r-process yields contribute to abundances as primary processes, and are difficult to 
separate. 
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2.5.6. P-Isotopes 

The z/p-process occurs in CCSNe environments where high neutrino fluxes create 
proton rich ejects via weak reactions with neutrons and protons (Frohhch et al. 2006). 
Additional neutrino interactions with the left-over protons after a-rich freeze-out produces 
a neutron abundance, which then undergo (n, p) strong reactions. The resulting protons 
can then be captured allowing the synthesis of proton rich isotopes up to mass number of 
A ~ 100 (Frohlich et al. 2006; Martinez-Pinedo et al. 2006). Since the evolution of this 
process depends upon //-interactions with free nucleons, it is independent of the initial star's 
metallicity, and is considered a primary process. 

Proton-rich isotopes beyond mass number ^ ~ 100 are believed to be created by the 
7-process, where successive photo-disintegration events on pre-existing metals occur in the 
interior layers of supernovae (Woosley and Howard 1978; Rayet and Prantzos 1990). The 
target nuclei for such events are previously synthesized metals from the r- or s-processes. 
This process generally begins with (7, n) reactions moving the target nuclei to the proton 
rich side of the stability line, where the rate of (7, p) and (7, a) start to dominate. 
Whereas all j>-isotope abundances are relatively small and do not dominate their respective 
elemental abundances, they cannot be accounted for by s- or r-process production, and so 
both of the j)-isotope production processes discussed are required to provide a complete 
galactic-chemical history. 

Similar to the main and weak components of the s-process, we expect to have a 
transition region between the z/j>-process and 7-process. A calculation of z/^^-process yields 
for a I5M0 star was performed by Thielemann et al. (2010). Their results show sufficient 
production of Mo and Ru isotopes relative to their solar abundances, which is much needed 
due to well known deficiencies in the 7-process production of these isotopes. From their 
results, the z/j)-process yields are shown to decrease quite rapidly beyond mass number 
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A = 100. This suggests that the region of overlap is fairly abrupt, with only ^°^Pd beyond 
a peak at ^^Ru owing non-negligible fractions of their abundances to the 7-process. 

A more complete analysis perhaps using a grid of stellar masses, in order to determine 
the distribution in the transition region, would likely offer only minor corrections to the 
final elemental scalings. Instead, new free parameters from the stellar and nucleosynthesis 
modeling would be introduced, however. The complete isotopic decomposition for the solar 
abundance pattern is shown in Appendix 1, and can be used for future isotopic reference. 

3. Model Description 

Our simple model tries to describe a typical average galactic composition, where the 
initial state is a homogeneous BBN composition. The final state is taken as a homogeneous 
composition equivalent to the isotopic solar abundance pattern. Each isotope from our 
decomposition of the solar abundance pattern is scaled as a function of a chosen normalized 
dimensionless parameter ^. At ^=0 the galaxy is in the BBN composition, and at ,^=1 
the galaxy has the solar composition. The benefit of this choice of parametrization is that 
its range is congruent with the total metallicity (relative to solar), Z/Z0. A comparison 
between ^ and Z/Z0 is shown in Section 5. This choice of parameter motivates choosing 
functional forms for the scalings of each process by addressing the predicted dependence 
on total metallicity these processes would obey. Note that this comparison should not be 
interpreted as a physical condition imposed upon this parameter. Indeed the parameter ^ 
is not a physical quantity. Each process responsible for isotopic production (and depletion 
for the case of D and ^He) is assumed to enrich (deplete) our model galaxy under the 
instantaneous mixing approximation, to preserve homogeneity across the entire metallicity 
range considered. 
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The parametrization of our model is an attempt to give isotopic evolutions characterized 
by typical trends in our Galaxy. For consistency we do not directly use metallicity as 
the argument of the evolutions; indeed metallicity is an output of the model, and hence 
cannot be used as an input. In practice, normalized metallicity Z/Z0 tracks ^ very closely, 
see Fig. 10, and the deviation is typically less than 2%. Whereas our model is not as a 
traditional GCE model that uses scaled stellar yields, we believe our parameterization 
allows us to improve upon the standard of scaling solar abundances by a constant factor, 
for use in nucleosynthesis studies. 

We proceed by addressing the functional forms for each process described in Section 2. 

3.1. Massive Stars 

The results of the massive star simulation (Heger and Woosley 2010) is used to model an 
isotopic abundance pattern at a low metallicity of [Fe/H] = —3, and represents a fixed point 
at that metallicity for the massive star isotopic scalings of our model. We use the common 
relation: [X] = Log(X/X0). Note that the "iron metallicity" (relative to solar), [Fe/H], 
should be distinguished from the total metallicity, [Z] : The former is the conventional proxy 
for the latter used in GCE models. For each isotope, the contributions from this model 
were interpolated linearly between their solar values to their respective abundances given 
by the simulation, at a parameter value of log(.^) = —2.5. The reasoning behind this choice 
of log(.^) = —2.5 is as follows. At solar metallicity massive stars are responsible for roughly 
30% of the total iron, but nearly 100% of the alpha isotopes. So then [Fejmassivc = —0.5, 
and [a]massivc = 0. At lowest metallicities, alpha isotopes effectively comprise the aggregate 
of metals by mass, since secondary sources and Type la SNe do not provide enrichment 
until later. Hence \Z\ ^ [a] massive holds in low metallicity regimes and will trail [Fe/H] by a 
nearly constant 0.5 dex until Type la onset. We then use our model parameter ^ in place 
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of "LJTaq. Below this parameter value of log(.^) = —2.5, each contribution was interpolated 
linearly from the fixed point to zero (to agree with the abundance pattern given by BBN). 
In doing this we assume that the abundances at [Fe/H] = —3 are from the first generation 
of stars, and below this the metallicity is due to varying amounts of mixing of the BBN 
composition. This comparison between ^ and [Fe/H] is approximate, and relies on the 
normalized metallicity Z/Z© tracking ^ very closely, see Fig. 10. 

3.2. Type la SNe 

There are three constraints for choosing a parametrization for Type la isotopes as a 
function of ^. First, Type la contributions experience a delay before they begin to enrich 
the ISM. This is due to the time necessary for, in the accreting white dwarf model, an 
intermediate to low mass star to evolve through its main sequences and then accrete 
sufficient material to surpass the Chandrasekhar mass limit. In the meanwhile, massive 
stars continue to enrich the ISM, making the galaxy more metal rich. The metallicity at 
which Type la are able to begin contributing to ISM enrichment is typically constrained to 
the interval —2 < [Fe/H] < — 1 depending on environment, galaxy size, etc., but is usually 
favored towards the upper bound. Hence, Type la contributions should be negligible below 
some value of ^. Second, upon crossing the Type la onset value, contributions should 
gradually rise to their solar values. Third, at high ^ the scalings should behave essentially 
like a primary process. Note that this last constraint is an assumption of our model and 
could be true only for a flat star formation rate (SFR). In fact investigation into Type la 
progenitors by Mannucci et al. (2006) predicts that the ratio of CCSNe to Type la rates to 
be increasing with redshift, a conclusion reached by posing a two component delay time 
distribution (DTD) progenitor model. If indeed there exists a Type la progenitor that 
operates at low redshifts and experiences a longer delay between the formation of the WD 
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and the later Type la explosion (named 'tardy' by Mannucci et al. 2006), a larger number 
of CCSNe relative to number of CCSNe at the birth of the WD contaminate the ISM, 
thus increasing the metallicity at a faster rate than it was increasing at the birth of the 
WD. Hence once the Type la SNe do explode, their products would have a metallicity 
dependence different than linear, due to this increasing rate ratio of CCSNe/SNe Type la. 
We are careful to note that it is unclear exactly how a non flat SFR would impact our 
scaling for Type la SNe, since the above discussion relies on time, and a relation between 
our model parameter and time cannot be formed. We simply note that there may indeed 
be an effect. The analogy presented above with respect to fiat SFRs serves to construct a 
rough behavior of Type la abundances, and does not mean to attach the physical quantity 
of time to our model parameter. Furthermore, Type la onset may not occur at one unique 
metallicity, there may be a spread - in how far the different environments have evolved in 
terms of Z (and O as its main tracer) - when Type la sets in. Hence, the corresponding 
parameter value ^ for la onset is not unique for all constituents, but rather just a typical 
average that fits the different environments that we sample as as function of it. 

One specific form which satisfies all three constraints discussed is given, e.g., by a 
scaled and shifted hyperbolic tangent base function, 

XYiO = ^fo ■ ^ ■ [tanh(a ■ ^ - b) + tanh(&)]/[tanh(a - b) + tanh(&)] (5) 

The specific value for Type la onset is determined by fitting the free parameters, a and 
b, against available data. Note the hyperbolic tangent function is tempered with a linear 
factor of ^ to ensure the behavior is linear near solar, where otherwise the tanh(x) function 
would asymptote. Note further that whreas this function is phenomenologically motivated 
it is not unique. The arctan(x) and erf(x) also satisfy the above constraints, however, the 
erf(x) is more computationally expensive than the arctan(x) and tanh(x) functions. 
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3.3. Neutron Capture and P-Isotopes 



The main s-process {Is, hs, and "strong" component), weak s-process, r-process, 
z/p-process, and 7-process contributions are parametrized as power laws, 
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The free parameters p and h denotes the power of the ^-dependence for primary and 
secondary {hs) processes, respectively. The weak component of the s-process was given its 
own parameter w. The motivation for choosing a power law dependence on these parameters 
lies in the relations between abundance and metallicity for primary and secondary events. 
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namely that primary events produce abundances linear in metallicity, and secondary 
produce abundances quadratic in metallicity. The chosen model parameter ^ takes the 
place of metallicity (since metallicity is indeed an output of the model and for consistency 
cannot be used as an input), and the exponents are parametrized for fitting with data. The 
metallicity dependence for the 7-process is less clear. It is possible for the target nuclei 
for this process to be either secondary or primary in origin, and produced in a previous 
astrophysical environment via the s- or r-process. Hence photo-disintegration events can 
potentially be a tertiary process. All 7-process isotopes have low (<1-10%) contributions 
to their respective elemental abundances, hence we chose not to assign a separate free 
parameter for this process, since the fit to data would likely be poorly constrained. Instead 
a compromise was struck between the possible s- or r-process origins and the chosen 
exponent for the 7-process represents an equal primary and secondary seed; about half the 
metals which might be target nuclei for this process are created by the s-process, and the 
other half from the r-process. Note that the proposed 7-process abundances enhanced from 
weak s-process seeds in the WSEP process (discussed at the end of section 2.4.) are scaled 
the same as the weak s-process, not as 7-process yields. 

The choice for the "strong" s-process function is motivated phenomenologically. At 
low ([Fe/H] < —2.0) metallicities essentially all of the solar abundances for the "strong" 
s-process have been made. Hence the abundances should be close to constant between 
—2.0 < [Fe/H] < 0. Additionally, since AGB stars experience a time delay in their 
contributions, below [Fe/H] < —2.0 the abundances drop smoothly, hitting zero at some 
unknown metallicity. The free parameters c, d, and g are used to adjust the shape of the 
hyperbolic tangent function, and constrain the peak abundance achieved at low metallicity, 
the metallicity at which the abundances begin to drop to zero, and the rate at which the 
abundances drop to zero. 
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3.4. Hydrogen Burning, Classical Novae, //-Process, and Galactic Cosmic Ray 

Spallation 

The scaling of traditional GCR spallation abundances has the same functional form as 
the s-process scaling. Novae, primary GCR spallation, up-process, deuterium, and helium 
(from hydrogen burning) abundances were scaled the same as the r-process, but with an 
offset added to reflect BBN abundances at .^ = 0. The remaining isotope of hydrogen, ^H, 
was scaled according to, 

XfiO = ^?0 ■ [1-0 - e ■ ^0 - YiO - DiO], (13) 

where Y{C,) is the mass fraction of the helium isotopes, D{C,) is the mass fraction of 
deuterium, and Zq is the solar value of total metallicity, given as Z0=O.O153 (Lodders et 
al. 2009). Effectively, this scaling is simply a restatement of the sum of mass fractions, 

X + Y + Z = 1. 

4. Fitting Scaling Model to Observational Data 

The isotopic scaling functions were summed into elemental scaling functions for the 
purpose of fitting the free parameters. A variety of compilations of stellar abundance data 
were used. The Frebel (2010) low metallicity data set contains over 1000 stars from the 
Milky Way (MW) halo and dwarf galaxies. The dwarf abundances were removed from the 
set, since it is likely that dwarf galaxies exhibit a different GCE than the MW. In fact a 
future extension of the current model could be applied to dwarf galaxies to give insight 
into their GCE relative to the MW. Additionally, Frebel (2010) data from stars in binary 
systems was removed from the [Ba/Fe] and [Sr/Fe] data, since binary systems experience 
more enriched s-process abundances due to accretion. Binary stars were identified from the 
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online CHARA catalog http://www.chara.gsu.edu/~taylor/catalogpub/catalogpub.html 
(Taylor et al. 2003), and were removed from the data set. In addition to removing the 
known binaries identified by the CHARA catalog, stars that simultaneously satisfied the 
following criteria: [Fe/H] < —2.0, [Eu/Fe] > 0.5, and [Ba/Fe] > 0.0 were also removed, since 
they are likely also candidates for binary systems not covered by the CHARA catalog. The 
other compilation used was the Soubiran and Girard (2005) data set, which contains, in 
part, 725 stars with magnesium abundances. The Soubiran (2005) data is needed for high 
metallicity abundances which fill a paucity of Mg data in the range — 1 < [Fe/H] < offered 
by Frebel (2010). 

All data is given in units of [X/Fe] as a function of [Fe/H]. The [Fe/H] axis was split 
into 300 bins in the range —5 < [Fe/H] < 1, and the [X/Fe] axis was split into 1500 bins in 
the range of the elemental data. Each data point was then assigned a Gaussian distribution 
using nominal values for the errors in the observations, 

fi = exp {-0.5 ■ ([(xi - Xo)/aj2 + [(y. _ Y,)/a,f) ] , (14) 

where fi is the value of the distribution at bin x; on the [Fe/H] axis and bin y; on the 
[X/Fe] axis, xq and yo are the [Fe/H] and [X/Fe] values of each data point, and a^ and 
ay are the nominal errors for each xq and yo. A nominal error of 0.1 dex was assumed for 
[Fe/H] and [Mg/Fe], and a nominal error of 0.15 dex was assumed for [Ba/Fe], [Sr/Fe], 
[Eu/Fe], and [0/Fe]. The Gaussian contributions were summed for each bin, and assigned 
an average and standard deviation for each bin. These binned averages and standard 
deviations were used for sampling the parameter spaces. We note that all fitting is done in 
the logarithmic space as is common practice. One could argue that abundances and yields 
are "additive," hence the linear space should be used. This is left for future work. Indeed 
fitting in the linear space would be impractical for large dynamic ranges on the one hand; 
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for ratios and small ranges the logarithmic space should suffice. 



4.1. Type la Parameters 



For fitting the Type la parameters a, b, and the fraction of the solar ^^Fe abundance 
attributed to Type la SNe, /, [Mg/Fe] data was chosen due the large number of data points 
that exist for this element as well as it owing its abundance to both massive stars and Type 
la SNe. The parameter spaces for a, b, and / were chosen to minimize x? between the 
[Mg/Fe] scaling model and the averages of the binned data. The results are given in Fig. 3. 
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Fig. 3. — The resulting model for [Mg/Fe] found by parameter fitting. Blue x's: Frebel 
(2010) data points. Yellow asteriks: Soubiran (2005) data. The dark shadow background 
shows the errors in the data, depicted by a Gaussian distribution about the data points. 
The averages and standard deviations for the Gaussian contributions per bin are given by 
the central and exterior pink lines, and the model itself is shown as a green line. There are 
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a couple outlying data points on the graph whose Gaussian contributions can be seen at the 
upper limit of the [Mg/Fe] axis. 

As shown in Fig. 3, the magnesium scaling falls well within the standard deviation 
of the data, and tracks the average well. The determined best fit parameter values are a 
= 4.640, b = 2.400, and / = 0.708, with a resuhing Xr = 0.0397. This low x? refiects the 
large spread inherent in stellar abundances and the observational uncertainties that exist, 
and also the fact that the bins (and data) are not uncorrelated, as each star is "spread" 
out over many bins due to the assigned Gaussians. The drop of the curve to solar values 
from its peak at [Fe/H] ^ — 1 is caused by Type la onset, when Fe production begins to 
dominate over magnesium. 

The fiat scaling depicted in the figure below [Fe/H] < —4 describes magnesium and 
iron evolving together. This is a consequence of our model assuming massive stars are the 
sole and unique source of metals at low metallicity. Whereas it is likely that magnesium 
and iron co-evolve in this range, due to their shared primary (massive star) origin that 
dominates at low metallicities, it is unsubstantiated from any data that the correlation is 
as exact as the model forces it to be. Indeed, the model cannot predict abundances below 
[Fe/H] ^ —4 with any reliability, given the paucity of data at such low metallicities. It 
could be argued that "average" scaling of abundances is not even a well defined concept 
in this range, since individual astrophysical events can have such a large stochastic effect 
on the metallicity content. In this sense, our model assumes "scaling" is a consequence of 
varying amounts of mixing with BBN from the same stellar sources, and until sufficient 
stellar events can reliably produce an average, our model is not a statistically accurate 
description of this mixing. 
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4.2. R-Process and hs Parameters 

For determining the value for the r-process parameter p, the chosen data was [Eu/Fe]. 
Europium is an r-process peak element with two isotopes, ^^^Eu and ^^'^Eu, both of which 
have dominant (~ 85 % its solar value) contributions from the r-process. We determined 
the optimized value for p using data from the Frebel set (2010). The best fitting scaling for 
[Eu/Fe] is given in Fig. 4. 
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Fig. 4. — The resulting model for [Eu/Fe] found by parameter fitting. The full scaling is 
shown in green. Blue x's: Frebel (2010) data points. The rest of the figure follows the 
convention of Fig. 3. 

The best fit parameter value was found to be p = 0.939 with Xr = 0.042. The 
previously found best fit values for a, b, and / (from Section 4.1) were used for Fe. Note 
that a nominal value of h= 1.5 for the s-process parameter was assigned for the purpose of 
fitting the [Eu/Fe] model. The choice of this nominal value for h has a negligible impact on 
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the best fit value for p due to the small s-process component of Eu, and for comparison an 
optimization of the parameter space for p with an h value of 2 yields a best fit value of p 
= 0.936 with x^ = 0.044, a difference of 0.3 % in p, and a difference of 4.7 % in x^- 

For fitting the values for the heavy main s-process parameter h, the chosen data was 
[Ba/Fe]. This element has two s-only isotopes, ^'^'^Ba and ^'^^Ba, along with three isotopes 
with contributions from both the s- and r-processes, ^'^^Ba, ^^^Ba, and ^^^Ba, and a small 
elemental contribution (~0.2 % its solar value) from two 7-process isotopes, ^^°Ba, and 
^■^^Ba. In Fig. 5, we plot results for the scaling of barium with the best fit value for h. The 
previously found best fit values (from Section 4.1) for a, b, /, and p were used for Fe and 
the r-process contributions to Ba. 
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Fig. 5. — The resulting model for [Ba/Fe] found by parameter fitting. The full scaling is 
shown in green, with the heavy s-process, r-process, and 7-process components plotted in 
teal, brown, and red, respectively. Blue x's: Frebel (2010) data points. The rest of the figure 
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follows the convention of Fig. 3. 

The minimized value of Xr = 0.134 was obtained from the parameter value h = 1.504. 
The trend of our barium model at low (< —2.0) metallicities implies abundances that 
begins to track iron. At these low metallicities, iron production is dominated by primary 
massive star contributions, as is barium production dominated by the primary r-process, 
hence the [Ba/Fe] encounters a "floor" in its abundances, as is expected (e.g., Truran 1981). 

Above [Fe/H] = —1.808, secondary contributions from the s-process exceed r-process 
contributions, and begin to drive [Ba/Fe] upward to a local maximum, before Type la 
contributions take the ratio back down to solar. This metallicity value of [Fe/H] = —1.808 
where heavy s-process contributions equal r-process contributions is lower than the typical 
value of [Fe/H] ^ —1.5 (found by, e.g., Travaglio et al. 1999; Truran 1981), and is due to 
our main heavy s-process exponent being smaller than the theoretical value of 2. Note that 
at all metallicities, the 7-process contributions to elemental abundances are negligible. 

4.3. Weak S-Process and Is Parameter 

The parameters constraining the weak s-process and Is scalings are w and /. Ideally 
one would wish to use [Ga/Fe] or [Se/Fe] elemental data, both of which have signiflcant 
weak s-process isotopic contributions to their elemental abundances (^ 0.61 for Ga, ~ 0.21 
for Se), as well as also lying on the flrst main s-process peak. Unfortunately data across a 
sufficient metallicity range for these elements is sparse, and [Sr/Fe] data is used instead, 
which has ^ 0.09 of its elemental abundance due to the weak s-process. In addition to 
Frebel (2010) data, two additional sources were used that provide observations of higher 
metallicity stars (Taylor et al. 2003). Sr has two s-only isotopes along the weak s-process 
path, ®^Sr and ®^Sr, along with one mixed isotope of r- and s-process origin, *^Sr (although 
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the r-process component is negligible), and one z/p-process isotope, ^^Sr. Possible binary 
contamination of the data was removed according to the same prescription adopted for the 
[Eu/Fe] data. In Fig. 6 we plot results of the model for strontium with the best fit value for 
w and /. The previously found best fit values for a, 6, /, and h (from Section 4.1 and 4.2) 
were used for Fe and the main s-process contributions to Sr. 
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Fig. 6. — The resulting model for [Sr/Fe] found by parameter fitting. The full scaling is 
shown in green, with the light s-process and weak s-process components plotted in cyan and 
brown, respectively. Blue x's: Frebel (2010) data points. Yellow crosses: Mashonkina et al. 
(2001) data points. Black asteriks: Jehin et al. (1999) data points. The rest of the figure 
follows the convention of Fig. 3. 

The value of Xr = 0.056 was obtained from the parameter value w = 1.230 and / 
= 1.227 . At all metallicities, contributions from the light s-process exceed weak s-process 
contributions, and the up-piocess and r-process contributions are negligible. 
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4.4. "Strong" S-Process Parameter 



The final free parameter constrains tlie tliird s-process peak. Pb data was taken from 
Frebel (2010), and binary stars were removed. The remaining data set consists of four 
points only. Due to the paucity of data, our usual standard of optimizing the elemental 
scaling using a x? analysis poorly constrains [Pb/Fe]. Instead, the functional form given in 
Section 3.3 was used to fit the free parameters c, d, and g by hand to the four data points. 
The result is given in Fig. 7. 
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Fig. 7. — The resulting model for [Pb/Fe] found by parameter fitting by hand. The full 
scaling is shown in green, with the heavy s-process and r-process components plotted in 
cyan and brown, respectively. Red x's: Frebel (2010) data points. 

As shown in Fig. 7, [Pb/Fe] peaks at [Fe/H] ^ —2.5, consistent with the AGB 
simulations by Bisterzo et al. (2010). Below [Fe/H] = —2.5, [Pb/Fe] drops until the 
"strong" s-process component vanishes, and the elemental scaling is determined only by the 
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r-process. The free parameters were found to be c= —2 ■ 10^^^, d= 200, and g= —0.23. It 
should be noted that the few data points can only constrain the peak of the [Pb/Fe] fit, and 
the drop of the abundance at [Fe/H] ^ —3.6 is not motivated by the data, but is an artifact 
of our chosen function for "strong" s-process evolution. It is clear that [Pb/Fe] should 
indeed drop, as AGB stars do not produce isotopes at arbitrary low metallicities, but the 
exact nature of the drop to the r-process may not be well represented by our model. 

Furthermore, the "kink" at [Fe/H] = —3 results from the change in slope of the Fe 
scaling. As discussed in Section 3.1, massive star contributions are linearly interpolated 
between their solar values and the values at [Fe/H] = —3 given by the massive star 
simulation. Below [Fe/H] = —3, massive star contributions are sent linearly to zero (in 
linear space). This treatment changes the slope of [Fe] on either side of [Fe/H] = —3, 
which manifests as the observed "kink" in Fig. 7. In reality we would expect a broad peak 
rather than the shown narrower peak following a kink before descent to the r-process 
floor. This very rough treatment of the "strong" s-process results in larger uncertainties for 
the s-process contributions to Pb and Bi isotopes from our model at metallicities below 
[Fe/H] < —2.5, and further revision is left to future work. 

5. Results and Discussion 

The functional forms for all scalings are now fixed by adjusting model to optimally 
fitting the observational data. A comparison is given in Fig. 8. 
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Fig. 8. — The scaling functions of the model contributions relative to oxygen as functions of 
metallicity. 

Type la SNe (yellow long-dashed line) contributions are negligible (< 1 % the 
solar value) below the Type la onset identified by our model. After Type la onset at 
[Fe/H] ^ —1.2, the contributions climb smoothly to solar. It is concerning that the Type la 
scaling begins to flatten at low metallicities until its slope becomes less than the 7-process 
(green dot-dot-dot-dashed line). Whereas contributions are negligible at these metallicities, 
the tanh(x) function chosen for Type la fails to describe the accurate physical picture, as the 
slope ideally should increase sharply below onset to reflect Type la "turning on." We accept 
this behavior for the Type la scaling in part because contributions are already negligible 
and would offer insigniflcant corrections to the isotopic abundances if flxed. Furthermore, 
the elemental data ends at [Fe/H] ^ —4, and hence the only constraint on Type la scaling 
below this metallicity (where it flattens out) is the BBN boundary condition. 
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Both the massive star (brown solid line) and r-process (black dotted line) scalings 
show similar trends at all metallicities, due to their shared primary nature. Abundances 
for the 7-process and z/p-process are scaled the same as the heavy s-process and r-process, 
respectively. 

Both components of the s-process do not show the typical behavior of secondaries. The 
heavy component (red dot-dashed line) shows a higher drop off at lower metallicities than 
the primary processes, but with a slower exponent of 1.504 compared to the theoretical value 
of 2. The light s-process (teal solid line) and weak s-process (pink dashed line) scalings 
behave similarly at all metallicities, with an exponent intermediate between the heavy 
s-process and primary processes. The "strong" component displays supra-primary behavior 
above [Fe/H] ^ —2.5, as is expected from its high abundance at low metallicities. Below 
[Fe/H] ^ —2.5 it decreases at a rate larger than any other process, until its contributions 
become zero and the scaling is no longer plotted. In this range [Fe/H] < —2.5 it is 
unlikely that our model correctly describes the "strong" component, as the scaling here is 
unconstrained by data. 

The failure of our model to reproduce the theoretical secondary nature is not new, 
and this discrepancy has been previously observed in the data (Prantzos 2011). It may be 
possible to alleviate this discrepancy by noting that at low metallicities rotating massive 
stars may have an increased neutron exposure in two different ways. The first is enhanced 
nitrogen production from CNO burning, which then burns into neon to seed the neutron 
source (Pignatari et al. 2008; Frischknecht et al. 2010). The second is an earlier ^^Ne 
ignition due to higher core temperatures (Frischknecht et al. 2010). The effect of this larger 
neutron exposure is enhanced weak s-process production at low metallicities, driving an 
increase in the abundance and resulting in a parameter value that is closer to a primary 
rather than secondary process. 
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Explaining the discrepancy for the heavy and hght components is more challenging. It 
may be plausible however, that we are observing the net result of a mixed isotopic history. 
That is, one could imagine a composition with a given abundance of an isotope made from 
a primary process in some astrophysical environment, that is then processed in a different 
astrophysical environment at a later time. If another abundance of this isotope is then 
made from a secondary process in the new astrophysical environment, its history would 
change from primary to secondary. Our model offers no way to track this effect. Indeed 
the aggregate of these effects may be in fact what we observe, and hence may be why our 
model gives an averaged main component parameter that deviates from the theoretical 
value of 2.0. Of course, this only applies to elements whose isotopic contributions can be 
made from both primary and secondary sources. In Table 1 of the Appendix we summarize 
our optimized parameter values. 



The complete scaling model for all elements is shown in Fig. 9. 



^^^^^^S^^^^^^S^^^^^^S^^^^^^H^^^^^^H^^^^^^H 
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Fig. 9. — The complete elemental scaling of the model. The abundances are given relative 
to their solar values. 

The traditional method of scaling solar abundances for inputs into stellar simulations 
is equivalent to treating the Galaxy as though all isotopic production is primary. A plot 
representing this approximation would be similar to Fig. 9, but with all metals changing 
their relative abundances at the same [Fe/H], which would look like a type of "flag" 
pattern, with all colors (representing relative abundances) changing together. In contrast 
to this approximation, Fig. 9 shows the corrections offered by our model to the traditional 
approximation due to the inclusion of Type la onset and secondary processes. These 
corrections are identifiable as "fingers" which protrude in the horizontal axis, distorting the 
otherwise clean "fiag" pattern, and occurring at elements that either lie on the Fe-peak or 
that have strong secondary or other Type la contributions. 

The scalings that notably stand apart from the others are the light elements: H, He, 
and Li. These begin with a much higher relative abundance (compared to the metals) due to 
BBN, and so change less, relative to their solar values compared to the other elements that 
begin with zero BBN contributions. Additionally, Pb displays a relative higher abundance 
at low metallicities due to the "strong" s-process. The sharp drop of Pb at [Fe/H] ^ —3.8 is 
consistent with Fig. 7, and is a very approximate treatment. 

Finally, our choice for the model parameter ^ can now be compared to the normalized 
metallicity Z/Zq. The plot is given in Fig. 10. 
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Fig. 10. — Comparison of our model parameter with 'L/'Lq. The behavior is hnear to within 
2% for all values of log(,^) and Z/Zq. The dashed lines demarcate the origin, and display 
the unity line for comparison. 

As shown in Fig. 10, our choice for the model parameter indeed corresponds to 
Z/Zq to within 2% for all values. The deviation from true linear occurs close to solar 
metallicity, when Type la SNe contributions become important and produce a small, almost 
unidentifiable "bump" near [Z] ^ —0.2. This comparison verifies our motivation for choosing 
functions for primary and secondary processes as being proportional to a polynomial of ^. 



Conclusions 



A metallicity-dependent Galactic isotopic decomposition for all stable isotopes has 
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been constructed. The solar abundance pattern was decomposed into several astrophysical 
processes responsible for isotope synthesis. Parametric functions were chosen to scale 
the contributions from each astrophysical process to give isotopic abundances, with the 
solar abundance pattern and BBN used as boundary conditions. The isotopic scalings 
were summed into elemental scalings and compared with stellar data in the halo and disk 
to tune the fit parameters of the model. The final scalings provide a complete isotopic 
abundance pattern at any desired metallicity. The purpose of this work is to provide 
isotopic abundances that can be used as initial abundances for stellar models in future work, 
or other nucleosynthesis studies. Our model is a large improvement over approximating 
the input isotopic abundance pattern by simply scaling the solar abundances, and/or using 
solar isotopic ratios. 

The interpretation of our offered solar abundance pattern decomposition is 
approximate, as different assumptions (of varying reliability) operate for different isotopes. 
The decompositions of the light isotopes (below carbon) likely reflect the best current 
understanding, however the current understanding is admittedly an area of ongoing 
investigation. The decomposition from carbon up to the Fe peak from massive stars and 
Type la SNe are the result of scalings done on data from two models that preserve the 
isotopic ratios across the two models, which distorts the original abundance patterns taken 
from the simulations. Whereas this preserves the salient features of these two processes, it 
should be noted that the isotopes that are not CNO or Fe peak suffer a larger uncertainty 
in their decompositions due to this scaling. Furthermore we make no distinction between 
the several operative primary processes that fall under the category of "massive stars," 
which are dominated by CCSNe yields but also necessarily include yields from stellar winds 
of lower mass stars, //-process, and classical novae. 

The main s-process yields used show good agreement with the solar abundances. 
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however, the weak s-process computation only reproduces the necessary abundances for 
three of the six s-only isotopes along the weak s-process path. The method for computing 
the weak s-process should ideally rely upon stellar models, unfortunately the results often 
display over-productions of several isotopes above their solar values. Our more simple 
approach, whereas less robust, has the advantage of only two isotopic over-productions, and 
is calculated directly from the neutron capture cross sections and branching ratios. 

In addition to providing isotopic abundances for input into stellar simulations, our 
methodology can also be applied to model other systems such as dwarf galaxies, and 
tailoring the free parameters to relevant data sets will then yield approximate chemical 
abundances for such systems. Yet another application lies with stellar model fits. Currently 
published stellar evolution models only give average yields for stars over certain mass 
ranges. These averages must be fitted to available observational data sets, as done for 
example by Heger and Woosley (2010). The need for this fit is evident, since stellar models 
only address the stellar nucleosynthesis part of GCE and neglect several processes which will 
ultimately infiuence the subsequent ISM abundances over many stellar lifetimes, processes 
such as infall and ISM mixing. The age of the Galaxy (or timescale of GCE) is much 
greater than the typical stellar lifetimes that contribute significantly to ISM enrichment; 
hence this fitting is required to give stellar yields a more precise physical meaning in a 
GCE context. Thus in addition to providing the initial isotopic abundances for input into 
stellar simulations, the abundances can also be used to fit the resulting stellar yields from 
the model in a consistent fashion. 

This work has laid out a basic method of isotopic decomposition as a function of 
metallicity based on elemental observational data and underlying nucleosynthesis processes 
for a complicated environment like the Galaxy. Future work should have a more detailed 
look at specific and less complicated environments like, e.g., dwarf galaxies that have 
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different contributions and would hence allow one to constrain model parameters more 
uniquely (such as Type la onset contributions). Another challenge will be to relate and 
identify the different nucleosynthesis processes with the principle components found in 
observational work like that of Ting et al. (2012), and help to improve such principle 
component analysis based on physical nucleosynthesis processes. 

We would like to thank Yong-Zhong Qian for useful discussions on nucleosynthesis, 
Nikos Prantzos for providing data from his own GCE model, and Anna Frebel for helpful 
correspondence concerning uncertainties in stellar abundances. 
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Table 1 summarizes the optimized parameter values found by fitting to data. 

Table 2 shows the solar abundance pattern decomposition for all stable isotopes into 
the various astrophysical processes employed by the model: Big Bang Nucleosynthesis, 
zz-process/primary galactic cosmic ray spallation/novae yields (together in a single 
category), secondary galactic cosmic ray spallation, massive star yields (includes CCSNe, 
stellar winds, //-process, and r-process contributions from carbon through zinc). Type la 
SNe yields, main s-process (which include all of "strong". Is, and hs components), weak 
s-process, z/p-process, 7-process, and the r-process (from zinc through uranium). The 
solar abundances in column 2 are from Lodders et al. (2009) and are in units of mole 
fractions. The various astrophysical processes in columns 3-12 show the fraction of the 
solar abundance attributed to each process, and these fractions can be used to decompose 
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Table 1. Optimized Parameter Values 



Parameter Best-fit Value Description 



a 


4.640 


b 


2.400 


f 


0.708 


P 


0.939 


h 


1.504 


I 


1.227 


w 


1.230 


c 


-2.e-ll 


d 


200 


9 


-0.23 



Type la tanh Scaling Factor 

Type la tanh Shifting Factor 

Fraction of Solar ^^Fe from Type la 

Primary Process Exponent 

/is-process Exponent 

/s-process Exponent 

Weak s-process Exponent 

"Strong" tanh Coefficient 

"Strong" tanh Scaling Factor 

"Strong" tanh Shift Factor 
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any desired solar abundance pattern. The fraction values are rounded to three significant 
figures. Note that the remaining helium not made in BBN is from hydrogen burning, which 
not explicitly shown in the the table. 
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